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Nearest-neighbor Heisenberg antiferromagnet on a face-centered cubic lattice is studied by extensive Monte 
Carlo simulations in zero magnetic field. The parallel tempering algorithm is utilized, which allows to over- 
come a slow relaxation of the magnetic order parameter and fully equilibrate moderate size clusters with up to 
N ~ 7 x 10 3 spins. By collecting energy and order parameter histograms on clusters with up to N ~ 2 x 10 4 
sites we accurately locate the first-order transition point at T c — 0.4459(f) J. 
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Geometric frustration generally denotes inability of 
a magnetic system to find unique classical ground state. 
It arises due to a competition between exchange inter- 
actions for certain types of magnetic sublattices. The 
wellknown examples are triangular and kagome lattices 
in two dimensions and pyrochlore and face-centered cu- 
bic (fee) lattices in three dimensions. Intriguing be- 
havior of geometrically frustrated magnetic materials 
have attracted a lot of experimental and theoretical at- 
tention in the past decade pQ. Frustrated properties 
of Ising antiferromagnet on an fee lattice (figG) have 
been recognized a long-time ago |3J . The case of vector 
(Heisenberg) spins has been investigated to a lesser ex- 
tent. Experimental realizations of fee magnets include 
type-I antiferromagnet UO2 E] and type-II antiferro- 
magnet MnO U3- 

In the present work we investigate a nearest- 
neighbor Heisenberg antiferromagnet on an fee lattice 
described by the Hamiltonian 
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where Si is a classical vector spin of unit length. Ev- 
ery spin interacts with 12 nearest neighbors separated 
by (±0, ±0, 0), (0, ±a, ±a) and (±a,0,±a). Frustrated 
properties of the model become apparent if one cal- 
culates the mean-field transition temperature T C MF = 
i|min{J q }|. The Fourier transform of exchange inter- 
action is 

J q = 4j[cos(q x a) cos(q y a) + cos(q y a) cos(q z a) 

+ cos(ffea) cos(&.a)] . (2) 
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Fig.f : Face-centered cubic lattice as a network of edge- 
sharing tetrahedra 



The minimum is reached at Q = (7r/o)(l, q, 0) with ar- 
bitrary q and all equivalent wave-vectors in the cubic 
Brillouin zone. This set includes, in particular, the type- 
I antiferromagnetic structure with Qi = (n/a)(l, 0, 0) 
and the type-Ill ordering with Q3 = [n/a)(l, |, 0). The 
type-I (type-Ill) structure becomes the only absolute 
minimum of J q if a weak second-neighbor exchange is 
added to Eq. Q with ferromagnetic (antiferromagnetic) 
sign, see, e.g., [Hj. 

Degeneracy of the nearest-neighbor model Q is a 
consequence of decomposition of an fee lattice into a 
network of edge-sharing tetrahedra, such that every site 
is shared between 8 tetrahedra. The classical ground 
state constraint consists, then, in a requirement of zero 
total spin for every tetrahedron. This yields an infinite 
number of collinear and noncollinear state with differ- 
ent periodicity, all of them having the same classical 
energy. The harmonic spin-wave analysis shows that at 
low temperatures thermal fluctuations select collinear 
states by an order by disorder effect pj. Such a result 
can be most easily understood by a method suggested 
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in Ref . [7] , where the effect of short wave-length thermal 
fluctuations is shown to lead to an effective biquadratic 
exchange between neighboring spins. For the classical 
spins on an fee lattice thermal fluctuations generate the 
following low-temperature correction to the free-energy: 

^ = -§E(V s ;) 2 - (3) 

(ij) 

Such a biquadratic interaction favors collinear spin ar- 
rangement. Two examples of collinear ground states 
include type-I spin structure S, = ecos(Qi • rj) and 
type-Ill configuration Sj = ecos(Q 3 ■ r ; + jir). Ad- 
ditional collinear ground states are constructed from 
the above two configurations by selecting crystal planes 
(parallel to one of the cubic axes) with the Neel type of 
spin order and rotating all spins in such planes by 180°. 
All obtained collinear states have exactly the same free- 
energy in the harmonic approximation due to an extra 
gauge symmetry of a quadratic Hamiltonian [xj. Their 
degeneracy is lifted by anharmonicities in the spin-wave 
Hamiltonian — the problem, which to our knowledge has 
not been considered analytically. 

Finite temperature transition of a type-I fee antifer- 
romagnet has been studied by the renormalization group 
approach |U| • Absence of stable fixed points within 
the e-expansion suggests a first-order transition driven 
by thermal fluctuations. The above calculations have 
been performed for the case when the type-I structure 
corresponds to the absolute minimum, that is the case 
of the spin model with a significant ferromagnetic 
second-neighbor exchange. In the nearest-neighbor case 
the anomalous contribution of thermal fluctuations is 
further enhanced due to extra soft modes. Therefore, 
the conclusion about a first-order transition should re- 
main essentially unchanged. 

Numerical Monte Carlo (MC) simulations of the fi- 
nite temperature phase transition in a nearest neighbor 
fee antiferromagnet have been performed by a standard 
single spin-flip technique ^3 ^2 ^] . First-order na- 
ture of the transition at T c ~ 0.45 J was unambiguously 
established from the finite size scaling of the peak in 
the specific heat and has been further supported 
by the energy histograms collected at the transition 
point 12 . As for the type of magnetic ordering at 
low temperatures there is no consensus among different 
authors. The first works have suggested the collinear 
type-I antiferromagnetic structure ^3 E] , though no 
results for the magnetic structure factor have been pro- 
duced. In the subsequent more detailed study ^21 the 
low-temperature phase of an fee antiferromagnet was 
described as 'a collinear state with glassy behavior.' Ap- 



parent difficulty in simulation of an fee antiferromagnet 
at low temperatures comes from the above mentioned 
degeneracy between various collinear states, which cor- 
respond to different local minima of the free-energy 
functional. The collinear states transform into each 
other under rotation of all spins in one crystal plane. 
The local minima of the free-energy are, therefore, sep- 
arated by rather large entropy barriers ~ L 2 , L being 
a linear size of the system. A single spin-flip MC tech- 
nique cannot produce appreciable moves in the phase 
space between distinct collinear configurations. As a 
result, the spin system hardly relaxes to the absolute 
minimum and the magnetic structure factor does not 
exhibit good averaging in a reasonable computer simu- 
lation time. 

In the present work we shall apply the novel ex- 
change MC method ^Sj for simulation of a Heisenberg 
fee antiferromagnet. This modification of the standard 
MC technique, also called parallel tempering ^3], has 
been developed for spin glasses, which is an outstand- 
ing example of hardly relaxing spin systems with nu- 
merous local minima separated by macroscopic energy 
barriers. In the exchange MC technique, several repli- 
cas of the spin system are simulated in parallel at a 
preselected set of temperatures. After a few ordinary 
MC steps performed on each replica, adjacent in tem- 
perature replicas are exchanged with a specially chosen 
probability ^21- The ensemble of parallel tempering, 
thus, includes two Markov processes: stochastic motion 
in the multidimensional phase space of the spin model 
and a random walk along a one-dimensional array of 
replicas. The second auxiliary process helps to dramat- 
ically decrease the correlation time for the main stochas- 
tic motion by repeatedly heating a given replica to high 
temperatures, where it quickly looses memory about the 
low-temperature magnetic structure. The exchange MC 
technique allows to equilibrate moderate size systems 
and has gained popularity in simulation of spin glasses 
[T5] and frustrated Ising models ^Hl- This method has 
not been applied so far to geometrically frustrated vec- 
tor spin models. 

We have performed exchange MC simulations of the 
model Q on finite clusters with periodic boundary con- 
ditions and N = 4i 3 sites for L = 4, 6, 8, 10, and 12 
(N = 256-6912). The highest temperature for the ex- 
change MC ensemble was selected at T up = 0.75 J > T c , 
where the system equilibrates fast with a standard MC 
technique. The lowest temperature in our simulations 
was 0.1J. The intermediate temperatures have been de- 
termined empirically starting with a geometric progres- 
sion in such a way that the exchange rates for replicas 
swaps are roughly uniform and exceed 0.1 for the largest 
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Fig. 2: Temperature dependence of the structure factor 
corresponding to the type-I antiferromagnetic ordering 
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Fig. 3: Temperature dependence of the collinear (ne- 
matic) order parameter 



clusters. In particular, in the vicinity of the first-order 
transition point the temperature steps have to be de- 
creased in order to avoid a bottleneck for replicas drifts 
across T c . In total we have employed Nt = 72 replicas 
for most cluster sizes. The MC procedure for simula- 
tion of individual replicas is the multiflipping Metropo- 
lis method, see, e.g., JT], which includes several flipping 
attempts (to = 5 in our case) on every spin under the 
same local field before moving to a next spin. This 
method is especially suitable for an fee lattice with a 
large coordination number z = 12, when calculation of 
local fields takes a significant part of the CPU time. Af- 
ter sweeping 2 times over the whole cluster, the replica 
exchanges have been performed. The above procedure 
constitutes one exchange MC step. All replicas were set 
in random initial configurations and equilibrated over 
10 5 exchange MC steps. Measurements have been per- 
formed for additional 5 x 10 5 -4 x 10 6 exchange MC 
steps. Equilibration of various physical characteristics 
has been checked (i) by absence of time evolution over 
at least a half of measuring time, and (ii) by compar- 
ison to the results for a different set of initial replica 
configurations obtained by gradual annealing of a single 
replica from high temperatures. 

During simulations of the exchange MC ensemble 
the internal energy, the specific heat, the magnetic 
structure factor, and the collinear order parameter have 
been measured. Results for the internal energy and the 
specific heat are identical with the previous results ob- 
tained by the standard MC method |11M12| and are not 
discussed here. The magnetic structure factor is given 
by a square of the antiferromagnetic order parameter: 

where (. . .) denotes a thermal average. For q = Q3 the 
structure factor exhibits a rather weak T-dependence 



and scales as 1/N with increasing cluster size at all 
temperatures. This generally indicates an absence of 
the corresponding order in the thermodynamic limit at 
any T. 

Results for the structure factor of the type-I anti- 
ferromagnetic ordering are presented in fig|2] Contri- 
butions from three noncquivalent wave-vectors of the 
type-I structure are added together. Symbols are used 
to distinguish different curves, while the lines are drawn 
through the actual data for Nt — 72 replicas. For the 
largest cluster with N = 6912 spins only Nt = 39 repli- 
cas down to T = 0.3J have been equilibrated. The 
remarkable feature of the presented data is an inverse 
finite-size scaling at temperatures below T c ~ 0.45J: 
the order parameter increases with the system size. 
The equilibrium sublattice magnetization deduced from 
S aa (q) is still significantly smaller than 1 even at 
T = 0.1J. This is a consequence of thermally excited 
stacking faults, domain walls and other defects in an 
ideal type-I antiferromagnetic structure. Linear size of 
such defects coincides with the lattice size. Their con- 
centration, therefore, decreases with increasing system 
producing a significant increase of the order parame- 
ter. Unfortunately, the lattice sizes are still too small 
to observe an asymptotic thermodynamic behavior for 
S m (q), though they definitely point at the spin order- 
ing with the wave- vector of the type-I antiferromagnetic 
structure. The data for two large clusters exhibit a clear 
jump at the transition temperature indicating the first- 
order transition. 

Measurements of S aa (q) cannot distinguish between 
presence of three domains of a single-g type-I structure, 
which is always the case for finite clusters, and a non- 
collincar triple-q spin state. In order to study this aspect 
of the magnetic ordering in an fee antiferromagnet we 
define a collinear order parameter. This is a single-site 
characteristics, which is insensitive to periodicity (wave- 
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Fig. 4: Distribution function (energy histogram) for sev- 
eral lattice sizes at the first order transition T c = 
0.4459J. Inset shows energy histograms for the largest 
cluster L = 18 (jV = 23328) at T — 0.4458J, full line, 
and T = 0.446J, dashed line. 

vector) of the magnetic structure but describes instead 
breaking of a spin-rotational symmetry. The collincar 
order parameter is given by a traceless second-rank ten- 
sor: 

i 

It vanishes for a noncollinear triple-g structure and has 
nonzero value for a collinear state. In the present case 
P a/3 is a secondary order parameter: it couples linearly 
to a square of the primary antiferromagnetic parameter. 
For an XY antiferromagnet on a checkerboard lattice, 
P a/3 is, however, a primary order parameter and char- 
acterizes a spin-nematic state [J]. 

In MC simulations, a square of the order parameter 
(JSJ is measured, which is given by 

^ 2 oi = ^£<s?sf^>-i- (6) 

The corresponding results are shown in fig|31 Since P c 2 ol 
is proportional to the fourth power of an antiferromag- 
netic order parameter, it does not show an apprecia- 
ble jump at the first-order transition temperature. The 
data for the collinear order parameter exhibit a remark- 
able lack of finite size scaling. This indicates that spins 
over the whole lattice are predominantly parallel or an- 
tiparallel to a certain direction. At T = 0.1J the aligned 
component of spin is (S a ) « 0.78. Thus, the com- 
bination of the structure factor results (fig0) and the 
collinear order parameter data (figOJl points uniquely to 
the collinear type-I antiferromagnetic order in a Heisen- 
berg fee antiferromagnet. 

Finally, we present results for the energy and the 
order parameter histograms (distribution functions) at 
the transition point. Since an fee antiferromagnet ex- 
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Fig. 5: Histogram for the magnetic structure factor at 
the first order transition T c — 0.4459J 

hibits a weak first-order transition, we find it more ad- 
vantageous to perform histogram collection with a sin- 
gle replica instead of setting up an exchange MC en- 
semble. Histograms have been collected utilizing a hy- 
brid MC algorithm: 3 multiflipping Metropolis steps 
are combined with 11 overrelaxation (microcanonical) 
updates JTj. On average 2 x 10 6 configurations have 
been generated to build one histogram. Its accuracy 
has been checked by comparing the final distribution 
function to intermediate distributions obtained with 
proper rescaling from 10 6 and 5 x 10 5 configurations. 
The quality of the data shown in fig^ is significantly 
higher than in the previous work |12|. which allows us 
to locate more precisely the transition point. The first- 
order nature of the transition is clearly observed from 
a double-peak structure of the energy histogram, see 
fig 0| Positions of the two peaks do not change signif- 
icantly with growing cluster size, while probability for 
intermediate states rapidly drops with increasing L. At 
T c w 0.4459 J, the probability weights in the two peaks 
for the largest L = 18 cluster (N = 23328) differ by 
approximately 15%. To demonstrate a strong tempera- 
ture dependence of the relative weight of two peaks we 
present on the inset of figQ] the distribution functions 
at T = T c ± 0.0001J. The system clearly spends more 
time in the upper (lower) peak above (below) the tran- 
sition temperature. We conclude, therefore, that the 
first-order transition in a Heisenberg fee antiferromag- 
net takes place at T c = 0.4459 ± 0.0001J. 

The distribution function for a square of the anti- 
ferromagnetic order parameter is shown in figEl The 
largest cluster also develops a double-peak structure 
characteristic for the first-order transition. Smaller clus- 
ters have, however, significantly wider distributions for 
the magnetic structure factor, which is not surprising in 
view of many local minima for different collinear states. 
It is interesting to notice that the jump of S aa (Q) at the 
transition point determined from the histogram for the 
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L = 18 cluster (figEJ is close to the jump observed for 
the L = 12 cluster in temperature scan (fig|21), though 
the latter lattice does not exhibit any visible double- 
peak structure at T c . This indicates that apart from a 
close vicinity of the transition point the L = 12 cluster 
may be very close to the thermodynamic limit. 

In conclusion, we have performed Monte Carlo sim- 
ulations for a finite temperature transition in a nearest- 
neighbor Heisenberg fee antiferromagnet. The obtained 
results clearly demonstrate a first-order transition into 
a collincar typc-I antifcrromagnetic structure due to an 
'order by disorder' effect. Entropy mechanism for selec- 
tion of the magnetic ordering suggests an interesting se- 
quence of phase transitions for a weak antiferromagnetic 
second- neighbor exchange < J' <C T c . The higher 
temperature transition from a paramagnetic state is de- 
termined by thermal fluctuations and takes place into 
the type-I collinear antifcrromagnetic structure. At low 
enough temperatures T ~ J' the energy selection over- 
comes the entropy effect and a second transition from 
the type-I into the type-Ill structure occurs. 
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hospitality of the Condensed Matter Theory Institute of 
Brookhaven National Laboratory during the course of 
this work. 
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